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Electric transport through circular graphene quantum dots: Presence of disorder 
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The electronic states of an electrostatically confined cylindrical graphene quantum dot and the 
electric transport through this device are studied theoretically within the continuum Dirac-equation 
approximation and compared with numerical results obtained from a tight-binding lattice descrip- 
tion. A spectral gap, which may originate from strain effects, additional adsorbed atoms or substrate- 
induced sublattice-symmetry breaking, allows for bound and scattering states. As long as the diam- 
eter of the dot is much larger than the lattice constant, the results of the continuum and the lattice 
model are in very good agreement. We also investigate the influence of a sloping dot-potential step, 
of on-site disorder along the sample edges, of uncorrelated short-range disorder potentials in the 
bulk, and of random magnetic-fluxes that mimic ripple-disorder. The quantum dot's spectral and 
transport properties depend crucially on the specific type of disorder. In general, the peaks in the 
density of bound states are broadened but remain sharp only in the case of edge disorder. 

PACS numbers: 73.22.Pr, 72.80.Vp, 73.22.-f 
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I. INTRODUCTION 

Recent advances in the fabrication of single layer 
graphene structures have facilitated the realization 
and application of graphene nanoelectronic devices. 
Graphene nanoribbons (GNR)i with constant or varying 
widths, down to sub-lOnm, have already been prepared 
and operated as field effect transistors^ Also, quantum 
dots (QD), which were either plasma etched or carved 
out mechanically have been investigated recently^— The 
usual way to confine QDs electrostatically by external 
gates, however, is in general not practicable in pristine 
graphene because of the gap-less bandstructure. Also, 
the charge carriers behave like Dirac-fermions and thus 
the occurrence of true bound states is affected by the 
Klein-tunneling mechanism^! 9 - The latter can be seen ex- 
perimentally by measuring the current flowing through 
a sample with a tunable potential barrier across the 
graphene sheet. 10 ~ 13 

Therefore, physical or chemical effects that arc able 
to open a gap in the energy spectrum of single layer 
graphene are of vital importance for further electric ap- 
plications. Among the proposed mechanisms for the cre- 
ation of a spectral gap, are size quantization in armchair 
GNR, magnetic interactions^ between the edge states, 
as well as application of external electric potentials along 
the sample edges^ in zigzag GNR. Other proposed mech- 
anisms that are effective also in broad graphene sheets 
are strain-induced gap openings ) 16 ' 17 substrate-induced 
band gap formation,— and chemical effects of adsorbent 
atoms and molecules.— A gap opening attributed to a 
breaking of the sub-lattice symmetry and detected in epi- 
taxial graphene grown on a SiC substrate was shown to 
produce a gap of about 0.26 eV.— 

In the absence of a spectral gap it was theoretically 
shown that an electrostatically confined QD can accom- 
modate only quasi-bound states.— ~— At the Dirac point, 
i.e., at energy E = where valence and conduction band 
touch, the electronic transport through QDs of certain 



shapes was also considered^ In this special situation, 
sharp resonances in the two-terminal conductance were 
predicted. However, in the presence of a spectral gap 
around the Dirac point, which in principle can be created 
by one of the mechanisms mentioned above, true bound 
states were obtained] 25 ' 26 Additional energy gaps due to 
Landau quantization can also be induced in quantum dot 
physics by the application of a strong perpendicular mag- 
netic field^^ 

In this paper, we investigate the electric transport 
through a circular electrostatic potential in the presence 
of a spectral gap. In mesoscopic physics, a positive scat- 
tering potential is usually called a quantum anti-dot, and 
a quantum dot when the potential is negative. Owing to 
the chiral symmetry of graphene, our theory is valid in 
both cases, and for the sake of simplicity we always call 
it a quantum dot. We calculate analytically the trans- 
port cross-section derived from the continuum Dirac for- 
mulation and compare the results with the two-terminal 
conductance obtained numerically from a tight-binding 
(TB) lattice model. In doing so, we also study the limit 
of validity of the continuum Dirac fermion description 
that is considered to be a good approximation of the low 
energy physics in the vicinity of the Dirac point. 

Real graphene samples are hardly perfect due to the 
presence of disorder induced by adsorben t 28 ' 29 atoms and 
molecules, which mainly affect the unsaturated dangling 
bonds at the sample edges, or due to bulk defects and 
ripples i^2r— Therefore, we investigate the influence of 
such disorder effects on the spectral and transport prop- 
erties of an electrostatically confined graphene quantum 
dot. In particular, we find that the impact of one-dimen- 
sional uncorrelated random disorder potentials, which 
only disturb the a-sublatticc sites at one edge and the b- 
sublattice sites at the other, causes changes of the quan- 
tum dot properties that are different from the case of 
short-range or random magnetic flux bulk disorder. In 
addition, we analyze what happens to the bound states' 
energies when the boundary of the electrostatic potential 



confining the QD varies over a length of several lattice 
constants. The effect of disorder on the single-particle 
states at the edges of graphene QDs has already been 
discussed recently^ 

The outline of the paper is as follows. In section [III 
we specify the continuum Dirac model and the TB lat- 
tice model for graphene, both include an electrostatic 
potential defining the QD. For a piecewise constant ra- 
dially symmetric potential, the Dirac equation is analyt- 
ically tractable. We discuss the character of the elec- 
tronic eigenstates occurring in various regions of an en- 
ergy versus potential diagram. In section IPTl] we obtain 
the energies of the bound states of the isolated QD and 
compare them with numerical results of the correspond- 
ing TB lattice model. In section IIV1 we add an environ- 
ment to the isolated QD and calculate both the scattering 
cross-section and, more specifically, the transport cross 
section and compare it to numerical calculations of the 
two-terminal conductance obtained with a transfer ma- 
trix method. Finally, section|V]is devoted to a numerical 
study of various types of disorder and its implications on 
the density of states in the vicinity of the Dirac point. 



II. GRAPHENE QUANTUM DOT WITH 
RADIALLY SYMMETRIC POTENTIAL 

A. Continuum model 

In the following, we recall the basic notions of the Dirac 
description of graphene. It stems from the low-energy 
expansion of the TB Hamiltonian around the K and K' 
points in the first Brillouin zone where the conduction 
and valence band touch. The wavefunction is then a 
four-dimensional spinor ( £+ £_ rj + r\_ ) , £ denotes the 
wavefunction in the valley K, and r\ that in the valley 
K' = — K and the indices +, — denote the sublattice (the 
unit cell contains two points). There is no coupling be- 
tween the valleys and thus we reduce the Dirac equation 
to a 2 x 2 matrix form for the valley K 



E-A-V(r) 

id x + d y 



id x - 

E + A 



8y 

-V(r) 






0. (1) 



The states r\± in the other valley K' are obtained by the 
transformation x — > — x and they follow a similar treat- 
ment as £±. In Eq. ([!]), E is the energy and V(r) the 
radially symmetric potential. Also, we introduced a con- 
stant mass term A that will account for a gap 2A in the 
energy spectrum. The gap is assumed to be substrate- 
induced as found in epitaxial graphene^ Above we have 
set Hvf = 1, however, in physical units, the Fermi veloc- 
ity vf = 3ta/(2h) is obtained from the TB model where 
a = 1.42 A is the carbon-carbon distance, and t w 2.7eV 
is the hopping integral. 

Since Eq. (JTJ) is rotationally invariant, we decompose 
the wavefunction £ into components with fixed angular 



TABLE I: Choice of the Bessel functions I m (x), J m (x), 
K m (x), N m (x) with H& 2 \x) = Jm(x) ±iN m (x) that de- 
scribe the radial part of the electronic wavefunctions inside 
and outside the quantum dot. 





inside r < R 


outside R < r 


I 


k' = in' 


k = in 




k' = a/A 2 - (E - U) 2 
£ J; ±(r) oc 7j±i/ 2 (kV) 


K = VA 2 - E 2 
i.j-±{r) cc K J±1/2 (K.r) 


II 


k' = ^(E-U) 2 - A 2 
£j;±(r) oc Jj ±1/2 (k'r) 


as in I 


III 


as in I 


k = VE 2 - A 2 
fr.±(r) oc H^\) 2 (kr) 


IV 


as in II 


as in III 



momentum 



&=(*) 



y e 'f(^D 



£j;±(r) 



(2) 



where J is half-integer. Then, for a piecewise con- 
stant V{r), the remaining radial equation is easily solved 
in terms of the Bessel functions I m (x), J m (x), K m (x), 
N m (x), and Hm (x) — J m {x) ±iN m (x). We consider 
here V(r) = UQ(R — r), where U denotes the potential 
in the QD and R its radius (measured in units of a) . The 
Bessel functions that describe the electronic states in- 
side and outside the dot are chosen according to whether 
the local energy is outside or inside the gap, i.e., if the 
wavevectors k' (for r < R) and k (for r > R) are purely 
real or imaginary, respectively, as summarized Table |TJ 

In the 'energy vs. potential' diagram, as shown in 
Fig. [TJ one distinguishes between four different behav- 
iors of the wavefunction. The diagram is symmetric un- 
der the transformations E — > —E and U — >• —U, and 




FIG. 1: (Color online) Energy domains for the QD states. I 
(white): no states. II (black, red): bound states. Ill (gray, 
blue): tunneling states. IV (light gray, green): scattering 
states. 



can be understood as follows. Everywhere in the dia- 
gram, the wavevectors inside and outside the dot satisfy 
k' 2 = (E - Uf - A 2 and k 2 = E 2 - A 2 , respectively. In 
domain I, one has \E — U\ < A and \E\ < A, and conse- 
quently k' 2 < and k 2 < 0. Therefore, both wavevectors 
kl and k are purely imaginary and the wavefunction de- 
cays inside and outside the QD. The domain II consists 
of all the E and U values that satisfy the \E — U\ > A 
and \E\ < A inequalities (k' 2 > and k 2 < 0), so that k' 
is real and k is purely imaginary. True bound states that 
oscillate for r < R and decay as exp (— r\/A 2 — E 2 ) when 
r — > oo, can exist therefore only in the energy domain II. 
For domain III one has \E — U\ < A and \E\ > A, i.e., 
k' 2 < and k 2 > 0. Here, the wavevector k' is imaginary 
and k is real, so that the wavefunction decays inside and 
oscillates outside the QD (tunneling regime). The do- 
main IV is given by \E — U\ > A and \E\ > A, i.e., 
k' 2 > and k 2 > 0. Therefore, both k' and k are real, 
and the wavefunction oscillates inside and outside the 
QD (scattering regime). 

In the domain I, the Bessel functions that satisfy 
the Dirac equation are A'±Ij±i/2( K ' r ) for r < i? and 
A±K J ± 1 / 2 (Kr) for r > R. Here, A'± and A± are the 
amplitudes of the wavefunction and n and k! are the 
wavevectors inside and outside the QD, respectively. The 
other two Bessel functions are excluded as possible solu- 
tions because they diverge for r = and r — > oo, leading 
to infinite densities at the origin and non-normalizablc 
wavefunctions. The amplitudes arc determined from 
the Dirac equation and from the matching conditions at 
r = R for the wavefunctions inside and outside the QD. 
The matching condition leads to an equation for A' ± and 
A± that has no nonzero real solutions for k and k! , and 
hence, there can be no states that decay inside and out- 
side the QD. In the following sections, wc discuss the 
domains II-IV. 



B. Lattice model 

The single band TB Hamiltonian that describes the 
non-interacting electrons in graphene in the presence of 
a radially symmetric potential reads 



H T b = J^(J7r + A r )4c r - t 



E 



etc 



(3) 



Here, c\ and c r are the fermionic particle creation and 
annihilation operators at the sites r of a hexagonal lat- 
tice with carbon-carbon distance a. The (r ^ r') denote 
nearest-neighbor sites and t = 2.7 eV is the hopping pa- 
rameter. In the following, we set the energy scale by 
putting t = 1. The potential is U r = U for sites inside 
the QD (r < R) and U r = outside. Also, we have in- 
troduced a lattice anisotropy A r for the two sublattices 
of graphene, A± = ±A. This anisotropy opens an en- 
ergy gap 2A in the spectrum and allows for the existence 
of true bound states. In the presence of disorder to be 



considered later, the above Hamiltonian is modified to 
account also for random on-site potentials or for random 
transfer terms between nearest neighbor sites. 



III. BOUND STATES OF THE ISOLATED QD 

A. Continuum model 

In this section, the energies of the bound states are 
calculated. To this end, we put V(r) = U<d(R — r) and 
consider the region II in Fig. [TJ The solution of the radial 
Dirac equation for r ^ R is given in terms of Bessel 
functions as 



£fc(r) 



A' ± J J±i (k'r), r<R 
A±Kj±i{nr), r>R. 



(4) 



The other two Bessel functions are excluded as solutions: 
for r —> 0, Nj±i(k'r) diverges and the density would be 
infinite at the origin, for large r, Ij±i(k,t) diverges and 
the wavefunction would not be normalizable. The ampli- 
tudes A- and A'_ are tied to A+ and A' + , respectively, 
by the equations of motion 




There are two matching conditions at r = R. Thus one 
obtains the following condition for the existence of a non- 
zero wavefunction (bound state) 



E-U + A Jj+^k'R) 
E-U - A Jj_i(k'R) 



A + E K j+ i(kR) 
A-E Kj_i(kR)' 



(6) 

The other valley yields a corresponding equation with 
(J + 1/2) replaced by —(J — 1/2). Since we have to 
consider all half- integer J's in Eq. ©, each bound state 
is doubly degenerate when one considers both valleys. 
In Fig. [21 wc show the evaluation of the energies of the 
bound states vs. the QD radius and vs. the value of the 
potential. Since everything is symmetric under E — > —E 
and U — > — U as mentioned above, we show only states 
for U > 0. For clarity, we present only the states with 
J = 1/2, 3/2, and 5/2. Here, we use A = 0.1 in TB 
units t = 1 (which corresponds to A = 0.27 cV). Note 
that when R — 25 a and U = 0.2, there are six bound 
states in the QD: three with J = 1/2, two with J = 3/2, 
and only one with J = 5/2, as shown in Fig. [2] When the 
QD radius or the strength of the confinement potential 
increases, more bound states can be accommodated in 
the QD. Our results are in excellent agreement with pre- 
vious calculations of bound state energies in a radially- 
symmetrio 2 ^ or a rectangular-shaped^ graphene QD. 
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FIG. 2: (Color online) The energies of the bound states in a 
graphene QD as function of the dot radius (left panel) and 
potential strength (right panel). We use A = 0.1 and show 
the results for J = 1/2 (solid line), J = 3/2 (dashed line), 
and J = 5/2 (dotted line) only. The intersections of the 
solutions with the vertical lines yield the bound states energies 
for R = 25a and 17 = 0.2. 
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FIG. 3: (Color online) Comparison between the eigenenergies 
of the TB Hamiltonian (circles) of the graphene lattice and 
the solutions of the Dirac equation (lines) as a function of the 
potential U. The parameters used are A = 0.1 and R = 25 a. 



B. Lattice model 

It is well known that the above results obtained within 
the commonly applied continuum Dirac-fermion approx- 
imation have a limited range of validity. To investi- 
gate this in more detail and to check possible devia- 
tions originating from neglecting lattice effects, we com- 
pute the eigenvalues of the lattice model via a numeri- 
cal diagonalization of the TB Hamiltonian ([3]). In our 
calculations, we study a graphene ribbon consisting of 
2N ac x N zz = 2 • 100 x 100 carbon atoms with periodic 
boundary conditions (N ac and iV zz are the number of 
armchair and zigzag lines, respectively). For the lattice 
anisotropy, we take A± = ±0.1 and calculate the eigen- 
values for several QD radii as a function of the potential 
strength. Changing iV ac or N zz does not alter the re- 
sulting eigenenergies E between —A and A as long as 
iV ac a « N zz a > R > a. 

Figure [3] exhibits the eigenvalues of the TB Hamilto- 
nian (|3]) vs. the potential U in the energy range \E\ < A. 
For comparison, we show also the eigenenergy solutions 
of the Dirac equation for the bound states as given by 
Eq. ©. We restrict the angular momentum to J < 7/2 
for the sake of clarity. There is an excellent agreement be- 
tween the two approaches for small angular momenta and 
small U. The angular momentum is not a good quantum 
number in the lattice model, nevertheless, we can identify 
the energies in both models for small angular momenta. 
If larger QD radii are considered, the good agreement 
between the TB and the Dirac equation becomes bet- 
ter even for higher values of J (results not shown here 
because there are too many eigenvalues in the figure, al- 
ready for R ~ 60 a). The reason for the slight deviations 
for small QD radii are the irregularities at the QD border 
appearing in the lattice model but which are not present 



in the continuum model. The corresponding lattice QD 
is not a perfect circle but has a chiseled edge. 



IV. TRANSPORT PROPERTIES OF A 
GRAPHENE QD 

A. Continuum model 

In this section, we want to calculate the transport 
properties through the circular QD. For this purpose, we 
need to introduce an environment to the isolated QD so 
that the exponentially decaying bound states are still fi- 
nite when reaching the outer region. We take a radially 
symmetric model with the following potential 



V(r) 



U, < r < R 
0, R<r <L 
V, L<r 



(7) 



The profile of the potential, which now consists of three 
regions, is schematically shown in Fig. 2] The radially 
symmetric choice makes the analytic determination of 
the S-matrix possible. Next, we explain our setup and 
sketch the derivation of the S-matrix. 

Compared to the isolated QD, there is a third region 
for the radial component, r > L. Here, the potential 
V is chosen such that the wavefunctions are oscillatory. 
Then, we can have incoming and outgoing waves in the 
environment and define a scattering cross-section. In the 
intermediate region, R < r < L, the wave functions de- 
cay, and inside the QD, < r < R, the energy can belong 
to bound states. This corresponds to the regime II of the 
isolated QD. 
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FIG. 4: (Color online) Potential landscape of a QD with ra- 
dius R and an additional environment for r > L. The vertical 
dotted arrows indicate \E — V(r)\ in each of the three regions. 
For R < r < L, the wavefunctions decay. Also, k", k' = in' , 
and k are the wavevectors inside the QD, outside, and in the 
environment region, respectively. 



Next, wc turn to the definition of the scattering matrix 
SM- ~— . For r — > oo, we write the wave function as 



1 



£t(x) = -=VF^±Ae' 



1 



y/ZTrkr 



VE-V±Af ± (<p) 



Akr 



(8) 



The first term is the incoming plane wave in x direc- 
tion, the second the outgoing spherical wave. Further- 
more, f±{<p) is the scattering amplitude, with f+{p>) = 
e lv f-((p), and ip is the scattering angle. The two com- 
ponents of £ individually solve the Dirac equation for 
graphene and are normalized such that the incoming 
(plane wave) particle current density is k, while the num- 
ber of particles leaving per unit time radially in the 
direction <p is dtp± (/*(^)e*>/_(^) + fl{<p)e~^ f+{<p)) . 
Thus, the differential scattering cross section is 

(9) 
We then expand in angular momenta with coefficients fj 

/±M = iE eW±iH ^ ( 10 ) 



v/2' 



j 



and define the scattering matrix Sj as Sj = 1 + fj to 
obtain 

Next we sketch the calculation of Sj. In our geome- 
try, the general solution of the Dirac equation for fixed 



angular momentum J is given by 



A±Hf ±h {kr) 



B±H { ?l x {kr), L<r 



■j±i 



A'i J, 



(k"r), 



B' ± / j± i(kV), R<r <L 
r <R. 



Here, the wave vectors are k = y/(E — V) 2 — A 2 , k' = 
y/A*-E 2 , and k" = ^J{E-U) 2 -A 2 . Using in © 

OC 

e ikrcos( v ) = J^ e i( ~ v+ ^ m J m {kr) , (12) 



m— — oo 



we can easily find the expression for A± and B± from the 
ansatz Eq. ((5J), i.e., the relation between A± and B± and 
5,/. Then, ,A± and B± are related to A'± and $±, and 
these in turn to A'± by the matching conditions. Thus, 
the elements of the scattering matrix Sj are determined. 
We skip the tedious but straight-forward (linear) algebra 
and just quote the result 



Sj = - 



detl?( 2 ) 
detXK 1 ) 



(13) 



and 
p(i,2) 
/ 






r(2,l) 



w' + Kj + (x' R ) w' + I J+ {x' R ) w>;j J+ (x")\ 
w'_Kj_(x' R ) -w'_Ij_{x R ) wlJj_(x") 



w + H)^>{x) w' + K J+ (x' L ) w' + I J+ (x' L ) 

{w-H^ix) w'_Kj_{x' L ) -xi/_Ij_{x' L ) 



/ 

(14) 



The appearance of the determinants is easily explained: 
they result from the solution of the 4 linear matching con- 
ditions (2 for £± at r = L and 2 more at r — R). Again, 

Hm 0*0, -Km 0*0 j Im(x), and J m (x) are the standard 
Bessel functions and we use the following abbreviations: 



w± 



:= y / \E-V±A\ , w± := Vl^-t/iAl 



:=V|-E±A|, J±:=J±-, 



kL, 



J R;L 



:= re'-R; ft'X 



:= k"R 



(15) 



To make connection with the two-terminal conduc- 
tance that can be measured in experiments, we calculate 
the transport cross section as^ 



a tr 



dtp (1 — cos if) 



da(ip) 
dip 



(16) 



For the scattering cross section from Eq. (|llj) we use the 
5j computed according to Eqs. (|13H15p . The results are 
shown in Fig. where we compare the bound state en- 
ergies of an isolated graphene QD of fixed size R and for 
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FIG. 5: (Color online) Transport cross-section of a graphene 
quantum dot (R/a = 30, A = ±0.1) with environment (upper 
panel). The position of the peaks correspond to the energies 
of the bound states of an isolated graphene QD with dot po- 
tential U = 0.2. They also agree with the results from the 
corresponding TB model (circles in the lower panel). The 
transport cross section (green dotted line) is shown in arbi- 
trary units with the 'background' value (red solid line) given 
by the empty dot U — 0. 



fixed potential U with the positions of the resonances in 
the transport cross section of the QD plus environment. 
The bound state energies and the positions of the reso- 
nances exactly coincide. When the energy of an incoming 
particle matches that of a bound state, the particle tun- 
nels through the QD and the cross section is small. The 
'background' cross-section is given by atr{U = 0), i.e., by 
the cross section of an 'empty' QD. 



B. Lattice model 

In this subsection, the two-terminal conductance g 
through a graphene quantum dot is calculated numer- 
ically within a lattice model and compared with the 
results of the transport cross-section of the continuum 
model. The two-terminal conductance^ is obtained from 
the transmission matrix T employing the transfer-matrix 
method^ according to 



9 = y TrTtT = ^(cosh 2 ( £j /2))-\ (17) 



where the Si parameterize the eigenvalues of T, and N 
is the number of open channels. The electron transmis- 
sion through a sample of width L y = 96 a and length 
L x = 192 a with two semi-infinite leads attached at x = 1 
and x = L x is calculated numerically. Periodic boundary 



FIG. 6: (Color online) The two-terminal conductance g (up- 
per panel) calculated within a TB model compared to the 
normalized transport cross section (lower panel) of a contin- 
uum Dirac model. 



conditions arc applied in the y direction. The sublatticc 
anisotropy, the dot potential and radius are A/i = ±0.1, 
U/t = 0.2, and R/a = 30, respectively. The outcome is 
shown in Fig.|BJ where the two-terminal conductance (up- 
per panel) and the relative transport cross-section (lower 
panel) are plotted as a function of energy. The latter is 
normalized according to <r tr (U = 0)/<Jtr(U) — 1 while the 
former is displayed on a log-scale in order to make the 
comparison easier. The reason is that in contrast to the 
transport cross-section of the infinite system considered 
in the continuum model, the two-terminal conductance, 
due to the finite sample length in the lattice model, does 
not drop to zero between the resonances that are associ- 
ated with the bound states. For the energetic positions 
of the resonances, an excellent agreement is obtained be- 
tween the two approaches. However, the agreement be- 
tween the two is good only as long as i? > a. When 
the QD radius is reduced, the positions of the trans- 
port peaks, calculated within the continuum model, differ 
from those of the lattice model due to the irregularities 
at the QD edges appearing in the latter approach. 

It is clear that the electrostatic potential confining the 
QD cannot be infinitely steep in real samples. Therefore, 
we also studied the effect of linearly sloping dot potential 
steps. We find that the electronic bound states of the 
QD shift linearly to smaller energies when the slope of 
the QD walls decreases. Therefore, it can happen that 
the bound states' energies fall below the lower edge of the 
gap, E = —A. Hence, in experiments the dot potentials 
should be fabricated as steeply as possible. 
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FIG. 7: (Color online) Influence of random magnetic flux 
disorder on the density of states of a circular QD with ra- 
dius 30 a and dot potential U/t = 0.2. The sublattice dif- 
ference is A/t = ±0.1 as before. With increasing disorder, 
f/(h/e) = 0.005 (red), 0.01 (green), 0.02 (blue), and 0.05 
(magenta), the peaks get broadened and start to overlap. 



FIG. 8: (Color online) The density of states p(E; W c d) within 
the gap arising from the sublattice difference A/t = ±0.1 in 
the presence of edge disorder of strengths Wed/t = 0.1 (red), 
0.2 (green), and 0.5 (blue). The bound states energies of the 
quantum dot (radius 30 a, potential strength 0.2 t) remain 
almost unaffected by the edge disorder. 



V. INFLUENCE OF DISORDER 

Real graphenc samples are usually subject to various 
disturbances and imperfections. It is known that the 
electronic properties essentially depend on the particu- 
lar type of disorder present. In the following, we study 
the influence of possible modifications that may affect 
the bound states of an electrostatically confined graphene 
quantum dot. The lattice model Hamiltonian of the un- 
perturbed system ([3]) is now modified to allow for di- 
agonal and non-diagonal disorder. The latter arises in 
graphenc when either ripples or random magnetic fields 
are present^ -33 ' 41 The impact of random on-site disor- 
der that can be due to adsorbed atoms and molecules 
attached to the bulk or to the dangling bonds along the 
edges of the graphene sheet, is described by diagonal ran- 
dom disorder potentials 



H 






(U T + A r + e r ) c\c T - t 



E 



(18) 



Again, t = 2.7 cV is the nearest-neighbor transfer en- 
ergy for graphenc, r denotes the sites of the carbon 
atoms, and (r ^ r') are pairs of those neighboring sites 
on a bricklayer lattice 4 ^ that arc connected by bonds. 
The sample width and the length are L y = 96 a and 
L x = 192 a, respectively. Periodic boundary conditions 
(BC) are applied in the x direction (length) and Dirich- 
lct BC in the y direction (width). As before, U r and 
A r define the quantum dot potential and the sub-lattice 
anisotropy, respectively. The on-site potentials e r de- 
scribe either the uncorrclated random bulk disorder or 
the one-dimensional edge-disorder. A box-probability 



distribution P(e r ) = 1/W is assumed for the random 
potentials — W/2 < e r < W/2, where W is a measure of 
the disorder strength. The case of ripples and corruga- 
tions is modeled by an equivalent random-magnetic-flux 
(RMF) disorder model leading to complex phase factors 
that are denned by ^2 $rr' = 2Tre<j> r /h, where the sum is 
taken along the bonds around a given plaquette situated 
at r with the random flux </> r piercing through it. With- 
out a QD, the chiral symmetry is conserved 41 - - — in this 
disorder model. The random fluxes are also drawn at ran- 
dom from a box distribution —f/2 < 4> r < f/2 with zero 
mean, and h/e is the magnetic flux quantum. The RMF- 
disordcr strength / varies between < f /(h/e) < 1. 

The result of random-flux disorder is shown in Fig. [7j 
where the density of states (DOS) of a graphene sheet 
with quantum dot and sublattice anisotropy is plotted 
for energies around the Dirac point. The sublattice dif- 
ference is taken to be A/t = ±0.1, the quantum dot 
potential is U/t = 0.2, and the QD radius is R/a = 30. 
The peaks of the bound states are still clearly visible 
in the case of small disorder strength f /(h/e) = 0.005. 
With increasing disorder, however, the peaks get broad- 
ened and disappear for f /(h/e) = 0.05. Therefore, the 
sharp bound energy states of a clean graphene QD can 
be destroyed in the presence of sufficiently strong ripple 
disorder. 

A completely different behavior is obtained in the case 
of one-dimensional edge-disorder only4^— This can be 
seen in Fig. [S] where again the DOS is shown for energies 
within the gap A/t = ±0.1 of a quantum dot with po- 
tential U/t = 0.2. With increasing edge-disorder strength 
from W^/t = 0.0 to W cd /t = 0.5, only the DOS of the 
peaks lying outside the gap are broadened and eventu- 
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FIG. 9: (Color online) The density of states p(E; W) show- 
ing the bound states of the quantum dot (radius 30 a, poten- 
tial strength 0.2 1) within the gap arising from the sublattice 
difference A/t = ±0.1 in the presence of on-site disorder of 
strengths W/t = 0.0 (red), 0.05 (green), 0.1 (blue), and 0.2 
(magenta). The sharp DOS-peaks broaden and finally disap- 
pear with increasing disorder strength. 



ally completely fill the energy regions between the sharp 
bound states that remain unaffected. Hence, electrostat- 
ically confined graphcnc QD do not deteriorate in the 
presence of edge-disorder and remain suited for experi- 
mental spectroscopic studies. This behavior is due to the 
exponential decay of the edge states which do not get 
mixed with the localized bound states of the QD even in 
the presence of additional edge-disorder potentials. 

When in addition to edge disorder the random on- 
site potentials are also allowed to affect the bulk sites, 
the DOS of the quantum dot's bound states get broad- 
ened with increasing disorder strength as can be seen in 
Fig. [9] As in the case of the random-flux disorder, but 
in stark contrast to the edge-disorder only situation, the 
sharp peaks broaden and finally disappear for disorder 
strength W/t > 0.2. A similar behavior is seen in the 
two-terminal conductance where, with increasing disor- 
der strength, the single sharp resonances are replaced by 
a cluster of resonances which essentially depend on the 
particular disorder realization. So, uncorrelated short- 
range random potentials of sufficient strength destroy 
the discreteness of the QD energies. In order to fab- 
ricate electrostatically defined graphene quantum dots, 



only samples that are clean, at least in the bulk area, 
can be used, whereas disorder at the edges is of minor 
importance. 



VI. CONCLUSIONS 

The electronic properties of circular graphene quantum 
dots, which can be created electrostatically in the pres- 
ence of a sublattice asymmetry, were investigated with 
the aid of two different models. First, we used a contin- 
uum model described by a Dirac-type equation that can 
be solved analytically, and second, we applied a tight- 
binding lattice model, which was evaluated numerically. 
The dependence of the electronic bound states and the 
transport cross-section on the electrostatic potential and 
size of the quantum dot were calculated. Both the spec- 
tra of bound quantum dot states and the electric trans- 
port through the QD show a very good agreement be- 
tween the two models when the radius of the QD is much 
larger than the carbon-carbon distance a. For smaller dot 
radii, the agreement is always reasonable for the low lying 
bound states, but gets worse at higher energies where in 
the continuum model the large angular momenta essen- 
tially contribute. In experiments, QDs of sizes 20-50 nm 
are fabricated^^ which means R/a ps 70—176. For such 
radii, the agreement between the two models is very good 
even for bound states with high values of the angular mo- 
mentum J. 

The results of the transport cross-section for a QD with 
attached environment calculated within the infinite con- 
tinuum model coincided with the numerically evaluated 
two-terminal conductance obtained for a finite graphene 
sample with attached semi-infinite leads described by a 
tight-binding lattice model. Within the TB model, we 
also studied the influence of a sloping dot potential, which 
caused the energy levels to shift to lower energies, and of 
bulk and edge disorder on the bound states of electro- 
statically defined graphcnc quantum dots. 

The presence of disorder severely influences the quan- 
tum dot's spectral and transport properties. Only in 
the case of one-dimensional uncorrelated random edge- 
disorder, do the peaks of the density of states correspond- 
ing to the bound states remain sharp. For on-site or 
random-magnetic-flux (ripple) bulk disorder, the peaks 
of the bound states in the DOS are broadened and fi- 
nally disappear with increasing disorder strength, and 
thus the quantum dot loses its experimentally observable 
characteristic spectral fingerprint. 
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